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Abstract. We investigate, by means of computer simulations, shapes of non- 
convex bodies that maximize resistance to their motion through a rarefied 
medium, considering that bodies are moving forward and at the same time 
slowly rotating. A two-dimensional geometric shape that confers to the body 
a resistance very close to the theoretical supremum value is obtained, improv- 
ing previous results. 



1. Introduction 

One area of investigation in contemporary mathematics is concerned with the 
search for shapes of bodies, within predefined classes, which permit the minimiza- 
tion or maximization of the resistance to which they are subjected when they move 
in rarefied media. The first problem of this nature goes back to the decade of the 
1680s, a time when Isaac Newton studied a problem of minimum resistance for a 
specific class of convex bodies, which moved in media of infinitesimal particles, rar- 
efied to such a degree that it was possible to discount any interaction between the 
particles, and in which the interaction of these with the bodies could be described 
as perfectly elastic collisions |2- More recently we have witnessed important de- 
velopments in this area with the broadening of study to new classes of bodies and 
to media with characteristics which arc less restrictive. However, almost without 
exception, the studies which have been published have given special attention to 
classes of convex bodies — the convexity of a body is a sufficient condition for the 
resistance to be solely a function of singular collisions. Even the various studies of 
classes of nonconvex bodies which have emerged, especially in the last decade, are 
based almost always on conditions that guarantee a single impact per particle — 
[H [3 m E] ■ Only very recently have there begun to emerge some studies supposing 
multiple reflections, as is the case of the work of Plakhov [10] • 

In the class of convex bodies, the problem is normally reduced to the minimiza- 
tion of Newton's functional — an analytical formula for the value of the resistance. 
But, in the context of nonconvex bodies, there is not any simple formula known 
for the calculation of the resistance. Even if it is extremely complex, in general, to 
deal analytically with problems of multiple collisions, for some specific problems of 
minimization the job has not been revealed to be particularly difficult, there even 
being some results already available [SI [9]. If, on the other hand, we consider the 
problem of maximization, then in this case the solution becomes trivial — for any 
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dimension, it is enough that the front part of the body is orthogonal to the direction 
of the movement. 

And what if the body exhibits, besides its translational movement, a sHght ro- 
tational movement? When we think of this kind of problem, we have in mind, 
for example, artificial satellites, of relatively low orbits, which do not possess any 
control system which could stabilize their orientation, or other devices in similar 
conditions. In this situation we imagine that, over its path, the device rotates 
slowly around itself. 

The problem of resistance minimization for rotating nonconvex two-dimensional 
bodies has already been studied in [TU]: it was shown that the reduction of re- 
sistance, as compared with the convex case, does not exceed 1.22%. In its turn, 
the problem of maximization of the average resistance of bodies in rotation is far 
from being trivial, in contrast with that which occurs when we deal with purely 
translational movement. This class of problems was, therefore, the object of study 
of the work carried out by the authors in [111 I12j : nonconvex shapes of bodies were 
investigated which would maximize the resistance that they would have to confront 
if they moved in rarefied media, and, simultaneously, exhibited a slight rotational 
movement. With the numerical study which was executed, various geometrical 
shapes were found which conferred on the bodies rather interesting values of resis- 
tance; but it was in later investigations, performed in the follow-up of this work, 
that the authors managed to arrive at the best of the results — a two-dimensional 
shape which confers on the body a resistance very near to its maximum theoretical 
limit. It is this latest result which now is presented here. 

2. Definition of the problem for the two-dimensional case 

Consider a disc in slow and uniform rotation, moving in a direction parallel to 
its plane. We will designate the disc of radius r by Cr and its boundary by dCr- 
We then remove small pieces of the disc along its perimeter, in an e-neighborhood 
of dCr, with £ € M+ of value arbitrarily small when compared with the value of r. 
We are thus left with a new body B defined by a subset of Cr and characterized 
by a certain roughness along all its perimeter. The essential question which we put 
is the following: up to what point can the resistance of a body B be augmented? 
More than getting to know the absolute value of this resistance, we are principally 
interested in learning what is the increase which can be obtained in relation to 
the smooth body (a perfectly circular contour, in this case), that is. learning the 
normalized value 

^ Resistance(g) 
^ ' ^ ' Resistance(a) ■ 

It is possible, from the beginning, to know some important reference values for the 
normalized resistance: R{Cr) ~ 1 and the value of the resistance R{B) will have 
to be found between 0.9878 ([I^) and 1.5. The value 1.5 will be hypothetically 
achieved if all the particles are reflected by the body with the velocity v+ (velocity 
with which the particles separate definitively from the body) opposite to the velocity 
of incidence v (velocity with which the particles strike the body for the first time) , 
v"*" = —V, a situation in which the maximum momentum is transmitted to the 
body. It is also possible for us to know the resistance value of some elementary 
bodies of the type B. This is the case, for example, of discs with the contour 
entirely formed by rectangular indentations which are arbitrarily small or with the 
shape of rectangular isosceles triangles. As was demonstrated in [T^], these bodies 
are associated with resistances, respectively, of i? = 1.25 and R = V2- 

Apart from being defined in the disc Cr, it is assumed that the body to be max- 
imized is a connected set B gM."^, with piecewise smooth boundary dB. Therefore, 
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let us consider a billiard in K.^ \ B. An infinitesimal particle moves freely, until, 
upon colliding with the body B, it suffers various reflections (one at least) at regular 
points of its boundary dB, ending up by resuming free movement which separates 
it definitively from the body. Denote by convS the convex hull of B. The particle 
intercepts the d{coTwB) contour twice: when it enters into the set conv_B and in 
the moment that it leaves. L = |9(convi?)| is considered the total length of the 
curve 9(conv_B) and the velocity of the particle is in the first and second moments 
of interception is represented by v and v"*", and x and the respective points 
where they occur. As well, the angles which the vectors — v and v+ make with the 
outer normal vector to the section of d{couyB) between the points x and x~^ are 
designated ip and ip'^ . They will be positive if they are defined in the anti-clockwise 
direction from the normal vector, and negative in the opposite case. With these 
directions, both tp as well as ip"*" take values in the interval [— 7r/2, tt/2]. 

Representing the cavities which characterize the contour of B by subsets f2i, 
ri2, . ■ ., which in their total make up the set conv_B \ B, the normalized resistance 
of the body B (equation ([T|)) takes the following form (cf. [Hj): 

being Lq = \d{convB) n dB\ the length of the convex part of the contour dB, 
Li = |5(convi3) nOi|, with i = 1, 2, . . ., the size of the opening of the cavity fi^, and 
R{Q.i) the resistance of the normalized cavity fij, in relation to a smooth segment 
of unitary size, with 
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(3) i2(f2i) = - / / {l + cos {ip^ {x, if) — (p)) cos ip dip dx. 

° J -1/2 J -IT /2 

The function t^"*" should be seen as the angle of departure of a particle which 
interacts with a cavity fi^ that has opening of unit size and is similar to fJ,;, with 
the similarity factor \/Li. 

From equation we understand that the resistance of B can be seen as a 
weighted mean i^iLi/L = 1) of the resistances of the individual cavities which 
characterize all its boundary (including resistance of the convex part of the bound- 
ary), multiplied by a factor which relates the perimeters of the bodies convB and 
Cr- Thus, maximizing the resistance of the B body amounts to maximizing the 
perimeter of convB (|9(conv_B)| < |9C,.|) and the individual resistances of the 
cavities fi^. 

Having found the optimal shape f2*, which maximizes the functional Q, the 
body of maximum resistance B will be that whose boundary is formed only by 
the concatenation of small cavities with this shape. We can therefore restrict our 
problem to the sub-class of bodies B which have their boundary integrally covered 
with equal cavities, and in doing so admit, without any loss of generality, that each 
cavity £1^ occupies the place of a circle arc of size e ^ r. As with Li = 2r sin(e/2r), 
the ratio between the perimeters takes the value 

\d{conyB)\ ^ sm{e/2r) ^ _ {e/rf 
^ ' \dCr\ e/2r ^ 24 ' 

or that is, given a body _B of a boundary formed by cavities similar to fi, from Q 
and ([4]), we conclude that the total resistance of the body will be equal to the 
resistance of the individual cavity f2, less a small fraction of this value, which can 
be neglected when £ <C r, 

(5) R(B)^Rin)-^-^Rin). 
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Thus, our research has as its objective the finding of cavity shapes fl which 
maximize the value of the functional ((31), whose limit we know to be found in the 
interval 

(6) 1 < sup^Rin) < 1.5, 

as is easily proven using ([2]): if is a smooth segment, (p~^{x,ip) = —tp and 
R{^1) = I /_^(y2 /-^/2 ^'^^ i'^'fi)) cos ip dip dx = 1; in the conditions of maximum 
resistance (p~^{x,ip) — ip, thus R{il) < | J^^^^^ J^^J^^^2cos(pd(pdx = 1.5. 

3. Numerical study of the problem 

In the class of problems which we are studying, only for some shapes of il which 
are very elementary is it possible to derive an analytical formula of their resis- 
tance as we saw in the rectangular and triangular shapes previously referred 
to. For somewhat more elaborate shapes, the analytical calculation becomes rapidly 
too complex, if not impossible, given the great difficulty in knowing the function 
ip^ : [—1/2, 1/2] X [— 7r/2,7r/2] [— 7r/2, 7r/2], which as we know, is intimately re- 
lated to the format of the cavity Q. Therefore, recourse to numerical computation 
emerges as the natural and inevitable approach in order to be able to investigate 
this class of problems. 

There have been developed various computational models which simulate the 
dynamics of billiard in the cavity. The algorithms of construction of these mod- 
els, as well as the those responsible for the numerical calculation of the associated 
resistance, were implemented using the programming language C, given the com- 
putational effort involved. The efficiency of the object code, generated by the 
compilers of C, allowed the numerical approximation of ^ to be made with a suf- 
ficiently elevated number of subdivisions of the intervals of integration — between 
some hundreds and various thousands (up to 5000). The results were, because of 
this, obtained with a precision which reached in some cases 10^®. This precision 
was controlled by observation of the difference between successive approximations 
of the resistance R which were obtained with the augmentation of the number of 
subdivisions. 

For the maximization for the resistance of the idealized models, there were used 
the global algorithms of optimization of the toolbox Genetic Algorithm and Direct 
Search" (version 2.0.1 (R2006a), documented in |T|), a collection of functions which 
extends the optimization capacities of the MATLAB numerical computation sys- 
tem. The option for Genetic and Direct search methods is essentially owed to the 
fact that these do not require any information about the gradient of the objective 
function nor about derivatives of a higher order — as the analytical form of the 
resistance function is in general unknown (given that it depends on ip^(x,(p)), this 
type of information, if it were necessary, would have to be obtained by numerical 
approximation, something which would greatly impede the optimization process. 
The MATLAB computation system (version 7.2 (R2006a)) was also chosen because 
it had functionalities which allowed it to be used for the objective function the 
subroutine compiled in C of resistance calculation, as well as the (p'^{x, (p) function 
invoked in itself. 

3.1. "Double Parabola": a two-dimensional shape which maximizes re- 
sistance. In the numerical study which the authors carried out in [HI [12], shapes 
of ri/ defined by continuous and piecewise differentiable / : [—1/2,1/2] — > M.^ 
functions were sought for: 

(7) ftf = {(x, y) : -1/2 < a- < 1/2, < y < /(x)} , 
with the interval [—1/2, 1/2] x {0} being the opening. 
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The search for the maximum resistance was begun in the class of continuous 
functions / with derivative /' piecewise constant, broadening later to the study of 
classes of functions with the second derivative /" piecewise constant. Not having 
been able with these shapes to exceed the value of resistance R = 1.44772, we 
decided, in this new study, to extend the search to shapes different from those 
considered in ((T]) . We studied shapes il^ defined by functions x of y of the following 
form: 

(8) = {(x, y):0<y<K -g{y) <x< g{y)} , 

where h > and g : [Q,h] — ^ M.^ is a continuous function with g(0) = 1/2 and 
5(/i) = 0. 

Similarly to the study which was carried out for the sets flf, in the search for 
shapes fl^ , the functions g were considered piecewise linear and piecewise quadratic. 
If in the classes of linear functions it was not possible to achieve a gain in resistance 
relative to the results obtained for the sets f2/, in the quadratic functions the 
results exceeded the highest expectations: there was found a shape of cavity fi^ 
which presented the resistance R = 1.4965, a value very close to its theoretical 
limit of 1.5. There were also carried out some tests with polynomial functions of 
higher order or described by specific conical sections, but, not having verified any 
additional gain in the maximization of resistance, it was decided not to report the 
respective results. There therefore follows the description of the best result which 
was obtained, encountered in the class of quadratic functions x = ±g{y). 

The value of resistance of the sets fi^ were studied, just as defined in ([8]), in the 
class of quadratic functions 

9h,0iy) = ay^ + f3y + 1/2, para 0<y<h, 

where h > and a = ^^^f^^^'^ (given that gh^pih) ~ 0). In the optimization of the 
curve, the two parameters of the configuration were made to vary: /i, the height of 
the dVl^ curve, and /3, in its slope in the origin {g'{Q)). In this class of functions the 
algorithms of optimization converge rapidly towards a very interesting result: the 
maximum resistance was reached with h = 1.4142 and /3 = 0.0000, and assumed 
the value R = 1.4965, that is, a value 49.65% above the resistance of the rectilinear 
segment. This result seems to us really interesting: 

(1) it represents a considerable gain in the value of the resistance, relative to 
the best result obtained earlier (in [HI [12]), which was situated 44.77% 
above the reference value; 

(2) the corresponding set Q,^ has a much more simple shape than that of set 
J7/ associated with the best earlier result, since it is formed by two arcs of 
symmetrical parabolas, while the earlier one was made up of fourteen of 
these arcs; 

(3) this new resistance value is very near to its maximal theoretical limit, which, 
as is known, is found 50% above the value of reference; 

(4) the optimal parameters appear to assume value which give to the set fl^ 
a configuration with very special characteristics, as in what follows will be 
understood. 

Note that the optimal parameters appear to approximate the exact values h — \f2 
and /3 = 0. The graphical representation of the function R{h, (3) through the level 
curves, figure [TJ are effectively in concordance with this possibility — note that the 
level curves appear perfectly centered on the {\/2, 0) coordinates; marked on the 
figure by "+" . 

Note also, in figure [2 the resistance graph R{h) for /S ~ 0, where it can equally 
be perceived that there is a surprising elevation of resistance when h \/2. Thus 
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1.41 1.412 1.414 1.416 1.418 1.42 



(a) 



(b) 



Figure 1. Level curves of the R{h,f3) function. 



the resistance of the ft^'^i^ cavity was numerically calculated with the exact values 
h ~ \f2 and /3 = 0, the result having confirmed the value 1.49650. 
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Figure 2. Resistance graphic R{h) for (3 = 0. 



The shape of the set D.^'^-f with h = \/2 and /? = is a particular case with 
which is associated special characteristics which could justify the elevated value of 
resistance presented. The two sections of the shape are similar arcs of two parabolas 
with the common horizontal axis and concavities turned one towards the other — 
see figure [H But the particularity of the configuration resides in the fact that the 
axis of the parabolas coincides with the line of entry of the cavity (axis of x), and 
that the focus of each one coincides with the vertex of the other. 

In order to be easily referred to, this shape of cavity (figure[3^) will be, from here 
on, named simply "Double Parabola" . Thus, in the context of this paper, the term 
"Double Parabola" should be always understood as the name of the cavity whose 
shape is described by two parabolas which, apart from being geometrically equal, 
find themselves "nested" in the particular position to which we have referred. 

Since the resistance of the Double Parabolas assumes a value which is very close 
to its theoretical limit, in a final attempt to achieve this limit, it was resolved to 
extend the study even further to other classes of functions g{y) which admit the 
Double Parabola as a particular case or which allow proximate configurations of this 
nearly optimal shape. In all these cases the best results were invariably obtained 
when the shape of the curves approximated the shape of the Double Parabola, 
without ever having overtaken the value R = 1.4965. 
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Figure 3. (quasi) Optimal 2D shape — the Double Parabola. 



4. Characterization of the reflections in the "Double Parabola" 

Eaeh one of the illustrations of figure S] shows, for the "Double Parabola" , a 
concrete trajectory, obtained with our computational model. 




(a) X = 0.45, (y9 = 75°. (b) a; = 0.45, = 55°. (c) a; = 0.45, ^ = 35°. 




(d) a; = 0.3, = 75°. (e) a; = 0.0, 93 = 35°. (f) a; = 0.48, v? = 5°. 

Figure 4. Example of trajectories obtained with the computa- 
tional model. 

It is comforting to verify that, with the exception of one trajectory, in all the 
others the particle emerges from the cavity with a velocity which is nearly opposite 
to that which was its entry velocity. This is the "symptom" which unequivocally 
characterizes a cavity of optimal performance. Even in the case of the trajectory 
of the illustration (f), the direction of the exit velocity appears not to vary greatly 
from that of entry. 

If we analyze the five first illustrations, we may verify that there exists some- 
thing in common in the behavior of the particle: in describing the trajectory, the 
particle is always subject to three reflections. This appears to be a determinant 
characteristic for the approximation of the two angles of entry and exit. Although 
this conviction is by nature essentially empirical, the results of the analytical study 
which the first author carried out in [S] are heading in the direction of confirming 
that one very significant part of the "benign" trajectories — those in which the 
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vectors velocity of entry and of exit are nearly parallel; we call them so because 
they represent positive contributions to the maximization of resistance — suppose 
exactly three reflections. In that study the author managed to demonstrate some 
important properties which help in the consolidating the numeric results presented 
for "Double Parabola" . In particular, it was shown that 

• There are no trajectories of fewer than 3 reflections; 



• For angles of entry outside the interval (— <poj Vo)j with ip^ = arctan 




~ 19.47°, all the trajectories are of 3 reflections; 
• In trajectories with 4 or more reflections, the angular difference is delimited 
by 2ipo: \(f - (f+\ < 2ipo. 

5. Other possible applications 

Besides maximizing Newtonian resistance, it is exciting to verify that the poten- 
tialities of the Double Parabola shape found by us could also reveal themselves to 
be very interesting in other areas of practical interest. 

Given the characteristics of reflection which the Double Parabola shape presents, 
we can rapidly conceive for it a natural propensity for being able to be used with 
success in the design of retroreflectors. Retroreflectors are devices which send light 
or other incident radiation back to the emitting source. Ideally the rctrorcflcctor 
should undertake this function independently of the angle of incidence, something 
which docs not happen with existing devices. 

As we have seen, the Double Parabola, although it docs not guarantee the perfect 
inversion of all the incident radiation, carries out this function with great success: 
it guarantees a good approximation of the directions of the incident and reflected 
flows for a significant part of the angles of incidence, and even for the rest it does 
not permit that the lag reaches elevated values. We predict, therefore, as very 
promising its possible utilization in the definition of new geometries for the optical 
elements which make up retroreflccting surfaces. 

Retroreflecting devices, although they may be used in a wide range of technolog- 
ical areas, as is the case for example of optical communications in open space, are 
used in a massive way in the automobile industry and in roadway signalling, and 
thus we can understand very easily their utility. In [S] is presented a exploratory 
study about a possible way to take advantage of our result. 

6. Conclusion 

In the continuation of the study carried out previously by the authors in [111 112] , 
with the work now presented it has been possible to obtain an original result which 
appears to us to have great scope: the algorithms of optimization converge for 
a geometrical shape very close to the ideal shape — the Double Parabola. This 
concerns a form of roughness which confers a nearly maximal resistance (very close 
to the theoretical majorant) to a disc which, not only travels in a translational 
movement but also rotates slowly around itself. In figure [5] one of these bodies is 
shown. 

Noting that the contour of the presented body is integrally formed by 42 cavities 
O with the shape of a Double Parabola, each one of which with a relative resistance 
of 1.49650, from jH) and (01) we conclude that R{B) = ^^^^^0^R{n) « 1.4951 is 
the total resistance of the body, a value 49.51% above the value of resistance of the 
corresponding disc of smooth contour (the smallest disc which includes the body). 
We know that if the body were formed by a sufficiently elevated number of these 
cavities, its resistance would even reach the value 1.4965, but the example presented 
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Figure 5. (quasi) Optimal 2D body. 

is sufficient in order for us to understand how close we are to the 50%, the known 
theoretical majorant. 

Taking into account that the majorant 1.5 is a theoretical result which only 
signifies the non-existence of 2D bodies (with boundary formed by identical cavities) 
that exceed this resistance value, the shape found by us could even embody an 
optimal solution. In confirming this hypothesis — that the value of resistance 
which characterizes the Double Parabola effectively embodies the maximum limit 
which it is possible to reach with real two-dimensional shapes — our result would 
earn redoubled importance. Although easy developments are not foreseeable, this is 
an important problem which remains open, waiting for future contributions which, 
if they do not outdo our result, will permit the reinforcement of our conjecture. 
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